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Abstract 

In quantum field theory, physicists routinely use "normal ordering" of operators, 
which just amounts to shuffling all creation operators to the left. Potentially confus- 
ing, then, is the occurrence in the literature of normal-ordered functions, sometimes 
called "Wick transforms" . We aim to introduce the reader to some basic results and 
ideas around this theme, if not the mathematical subtleties; our intended audience 
are instructors who want to add something to their quantum field theory course, 
or researchers who are interested but not specialists in mathematical physics. For 
rigorous proofs and careful mathematical discussions we only give references. 
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1 Introduction 



Normal ordering was introduced in quantum field theory by G.C. Wick in 1950, to avoid 
some infinities in the vacuum expectation values of field operators expressed in terms of 
creation and annihilation operators The simplest example of such infinities can be 
discussed based on non-relativistic quantum mechanics and the simple harmonic oscillator; 
an infinite number of harmonic oscillators make up a free quantum field. (The reader who 
needs a quick reminder of some basic quantum field-theoretical concepts may find comfort 
in appendix A). We will use the harmonic oscillator to exhibit the connection between 
Wick-ordered polynomials and the familiar Hermite polynomials. Then we turn to Wick 
transforms in the functional integral formalism of field theory, where we show that there 
is again a connection with Hermite polynomials. Several different approaches to Wick 
transforms that can be found in the literature are compared and we show why they are 
equivalent. In passing, we observe how the standard quantum field theory result known as 
"Wick's theorem" follows rather directly in this framework, from well-known properties 
of Hermite polynomials. Finally, we provide one brief example of how the Wick transform 
can be utilized in a physical application. 



2 Wick operator ordering 
2.1 Simple Harmonic Oscillator 

The Hamiltonian operator for the simple harmonic oscillator in non-relativistic quantum 
mechanics has the form 

where we have, as usual, hidden Planck's constant h, the mass m, and the angular fre- 
quency u in definitions of dimensionless momentum and position operators 



1 ^ ^ fn^ ^ 



\/hmui V ^ 

so that 

[P,Q] = -z. 

If one defines the creation and annihilation operators and a by 



at = -^iQ-^P), (1) 
a = -^(Q + tP), (2) 



2 



so that 

[a,at] = l, (3) 

one finds that 

H = 1{P^ + Q^) = ^[a^a + aa^) = a^a + l . (4) 

This means, as proved in introductory quantum mechanics books (see e.g. Sakurai 0), 
that the eigenvalues of the Hamiltonian operator come in the sequence 

En = n + ^ , n = 0, 1,2, . . . . 

In particular, the ground state energy (or zero-point energy), which is the lowest eigenvalue 
of the Hamiltonian, is non-zero: 

H\0) = l\0). 

This is in agreement with Heisenberg's uncertainty principle: It is the smallest energy 
value that saturates the uncertainty relation (again, see e.g. Sakurai for the explicit 
calculation). 

This zero-point energy has observable physical consequences; as an illustration, it is 
possible to measure zero-point motion (which leads to the zero-point energy) of atoms in 
a crystal by studying dispersion of light in the crystal. Classical theory predicts that any 
oscillations of the atoms in the crystal, and therefore also dispersion effects, cease to exist 
when the temperature is lowered towards absolute zero. However, experiments demon- 
strate that dispersion of light reaches a finite, non-zero value at very low temperature. 

In quantum field theory, a free scalar field can be viewed as an infinite collection of 
harmonic oscillators, as described in appendix A. If we proceed as before, each oscillator 
will give a contribution to the zero-point energy, resulting in an infinite energy, which 
seems like it could be a problem. 

One way to remedy the situation is to define the ground state as a state of zero 
energy. We can achieve this by redefining the Hamiltonian: We subtract the contribution 
of the ground state and define a so-called Wick- ordered (or normal-ordered) Hamiltonian, 
denoted by putting a colon on each side, by 

: H: = : ^ (a^a + aa^) : =^ ^ (a^a + aa^) — |(0|a^a + aa^|0) = a^a. (5) 

Hence, in this example the definition of Wick ordering can be thought of as a redefinition 
of the ground state of the harmonic oscillator. 

On the other hand, we also see that in the last equality in (||) all creation operators 
end up on the left. A general prescription for Wick ordering in quantum field theory in 
a creation/annihilation operator formalism is then: "Permute all the and a, treating 
them as if they commute, so that in the end all are to the left of all a." The resulting 
expression is, of course, the same: 

: H: = a^a . 
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2.2 Wick ordering and Her mite polynomials 

The first connection between Wick ordering and Hermite polynomials arises if we study 
powers of the (dimensionless) position operator Q. Physically, in the harmonic oscillator, 
the eigenvalue of gives the variance (squared standard deviation) of the oscillator from 
rest. We have Q = {a^ + a)/V2, but to avoid cluttering the equations with factors of V2, 
we will study powers of just {a'^ + a). 

(a^ + a)^ = a^^ + a^a + aa^ + 

= a)'^ + 2a) a + + [a, a^] 
= : (a^ + a)^: + [a, a^] 

= -.{a^ + af-. + l by (5). 

Arranging terms in a similar way for higher powers of {a^ + a) we find 

(a'f + a)^ = : (a"^ + a)^ : + 3(a^ + a) . 

(flt + a)^ = : (a"f + a)^ : +6 : (a"f + a)2 : +3 . 

We can summarize the results as follows, with the notation + a = g, 

g3 = ■.q^:+ 3: q: , 

g4 = .g4. ^ Q:q^:+ 3. 

Since we can recursively replace normal-ordered terms on the right by expressions on the 
left which are not normal-ordered (e.g. : g^: can be replaced by — 1), we can also invert 
these relations: 

:q^: = q^-l = He2(g) , 

: g3 : = g3 _ 3^ = He3(g) , 

: g4 . = g4 _ g^2 ^ 3 ^ }iei(q) , 

where the polynomials He„(g) are a scaled version of the more familiar form of the Hermite 
polynomials H„: 

He,(x) = 2-"/2H„(x/y2). (6) 

Confusingly, in some mathematical physics literature, the He^ are often just called H„. 
Some of the many useful properties are collected in appendix B for easy reference (a more 
complete collection is given in e.g. Abramowitz & Stegun |^). 

Because of this relation between operator Wick-ordering and Hermite polynomials, the 
mathematical physics literature sometimes defines 'Wick ordering" in terms of Hermite 
polynomials: 

■.q'^-.'^Reniq). (7) 
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Although q is an operator composed of noncommuting operators a and a^, this alterna- 
tive definition naturally generalizes to Wick-ordering of functions. As promised, we will 
explore this idea in the next section. 

One reason that the connection to Hermite polynomials is not mentioned in standard 
quantum field theory literature is the fact (which was also Wick's motivation) that the 
normal-ordered part is precisely the part that will vanish when one takes the vacuum 
expectation value. Indeed, the traditional way to define normal ordering, the one given 
at the end of section 2.1 ("put to the left of a"), yields for powers of q simply 



which vanishes for any nonzero power n when applied to the vacuum state |0). 

In other words, since one knows that normal ordered terms vanish upon taking the 
vacuum expectation value, one may not be interested in their precise form. 

However, when the expectation value is not taken in the vacuum (for example, in a 
particle-scattering experiment), this part does of course not vanish in general, and there 
are in fact many instances where the actual normal-ordered expression itself is the one of 
interest . 

3 Functional integrals and the Wick transform 

Most modern courses on quantum field theory discuss functional integrals (sometimes 
called "path" integrals, however only in nonrelativistic quantum mechanics does one re- 
ally integrate over paths). In a functional- integral setting, the counterpart of the Wick 
ordering in the operator formalism is the Wick transform. This transform applies to func- 
tions and functional. It can, like its quantum-mechanics counterpart eq. (|^), be defined 
by means of Hermite polynomials. But first, let us briefiy skip ahead and explain why 
such a transform will prove to be useful. 

3.1 Integration over products of fields 

In the functional formalism, physical quantities like scattering cross sections and decay 
constants are computed by integrating over some polynomial in the fields and their deriva- 
tives. The algorithmic craft of such computations is described in textbooks such as Ryder 
|§] and Peskin & Schroeder P]. Even though there are examples of physical effects that 
can be studied with functional-integral methods but not with ordinary canonical quanti- 




n 



i=l 
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zationQ, within the scope of this paper we can only give examples of some things that can 
be derived more quickly or transparently using functional integrals. 

In fact, we will also be concerned with some rather basic questions that are usually 
glossed over in introductory treatments: what does the functional integral itself really 
mean? While a complete answer is not even known, and certainly beyond the scope of 
this short article, we intend to give some flavor of the first steps towards addressing this 
question and how the Wick transform has been put to work in this regard. 

First, a restriction: the polynomials considered in this section are polynomials of Eu- 
clidean fields (fields defined on four-dimensional Euclidean space R^). Similar formalisms 
exist for Minkowski fields (fields defined on spacetime) with minor changes in the equa- 
tions (see e.g. the aforementioned textbook Unfortunately, functional integrals over 
Minkowski fields are less mathematically developed than integrals over Euclidean fields, 
therefore we shall restrict attention to Wick transforms of functions and functionals of 
Euclidean fields — primarily polynomials and exponentials. The Wick transform, like the 
Hermite polynomials, has orthogonality properties that turn out to be useful in quantum 
field theory, as we shall see. First, we have to introduce a few mathematical concepts. 

3.1.1 Gaussian measures 

Here, our aim is to fix the notation, and to briefiy remind the reader how to integrate 
over Euclidean fields, without going into too much detail. The standard mathematical 
framework to perform such integration is the theory of Gaussian measures in Euclidean 
field theory, for which details can be found in the mathematical physics literature, such 
as Glimm & Jaffe [^] and Janson 

As a first try, one would define a field in the functional integral as a function on R^. 
Fields in the functional integral, however, may seem like functions at first glance, but can 
produce divergences that cannot (for instance) be multiplied in the way that functions 
can. A more useful way to regard a quantum field in the functional integral formalism is 
as a distribution $ acting on a space of test functions /: 

$(/) = ($,/), (8) 

where the bracket ( , ) denotes duality, i.e. $ is such that it yields a number when 
applied to a smooth test function /. In many situations the distribution $ is equivalent 

^In particular, contributions to the S'-niatrix that have an essential singularity at zero coupling con- 
stant cannot be found by standard perturbative expansion around zero coupling, yet these "nonpertur- 
bative" contributions can be studied using functional integrals. Rates for decays that would be strictly 
forbidden without these effects can be computed, see for example the book by Ryder Ch. 10.5]. 
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to a function 0, which means this number is the ordinary integral 

($,/)= [ d'x<j>{x)f{x). (9) 



A famihar example of a distribution is the one- dimensional Dirac distribution $ = 5, for 
which we have 

Hf) = {6J) = f{0). 

Now, just like a function, $ in general belongs to an infinite-dimensional space. To 
be able to integrate over this space (not to be confused with the integral in equation (H), 
which is an ordinary integral over spacetime) we need a measure, some generalization of 
the familiar dx in the ordinary integral above. Here, a useful generalization will actually 
depend on the Green's function, called covariance in this context and denoted C. In 
general, the covariance is a positive, continuous, non-degenerate bilinear form on the space 
of test functions. In the following, we will often encounter the covariance at coincident 
test functions, here denoted C(/, /). 

To get to the point, a Gaussian measure dfic is defined by its covariance C as 

j^df^cm exp(-^(<|.,/)) =exp(-lC(/,/)) (10) 

over a space Y of distributions $. For comparison, the usual Gaussian measure on M'^ is 
defined by | 



(27r)'^/2 J^, (detD)V2 
where 

{x',x)M.d = x'^x^ is the duality in W^, x G R'^, x' G Rd, 

Q{x) is a quadratic form on W^, Q{x) = D^i^x^^x" = {Dx,x)^d, 
W{x') is a quadratic form on Rd, W{x') = x'^C^^x^ = {x', Cx')^d, 

such that 

DC = CD = 1. (12) 

^ To avoid dimension-dependent numerical terms (powers of tt, powers of 2) in the definition ( |Il| ) of 
the Gaussian measure one can, alternatively, define it by 



/ 



(deti?)i/2 



g — 7rQ(x)g — 27ri{x',x) _ ^ — TrW{x') 



In fact, this can even be convenient in the simplest Fourier transforms for those who forget where the 2tt 
goes: f{p) = J dx e^'^'^'^P^ f (x) yields an inverse f{x) = J dp e'^^^^^ f{p), without forefactor. 
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In the more familiar M case (|TT]), it is the combination of standard measure and kinetic 
term that corresponds to the measure dfic we introduced above: 

e"2<3(^) is analogous to d/ici^) , (13) 



(detD)V2 



but going back from dfic to an explicit separation as on the left will not turn out to be 
necessary for our discussion. In fact, by defining the measure dfic through eq. ([T0|), we 
have not even specified what such a separation would mean. 

With the above expressions in mind, the covariance at incident points is expressed as 
the following integral, obtained by expanding eq. (p!0|): 

c{fj) = J dficm {^jf. (14) 

In fact, the integral on the left-hand side of eq. (0) is the generating function of the 
Gaussian measure; let us denote this integral by Z{f). This means that by successive 
expansion of eq. (|10|), the n-th moment of the Gaussian measure can be compactly written 
as 

I rf/ic($) /r = 

in-m CifJY/^ n even 

(15) 

n odd, 

where nil = n{n — 2)(n — 4) ■ ■ ■ is the semifactorial. For convenience, we introduce the 
following notation for the average with respect to the Gaussian measure ^c'- 

{Fmm,a--= J dficmrmf)]. (16) 

Note the difference between the brackets ( used for average and the brackets ( , ) 
used for duality. Armed with this set of definitions, we can define a Wick transform of 
functionals of fields. 



A=0 



3.1.2 Wick transforms, definitions 

The goal here is to provide some idea of how to address the difficult mathematical problem 
of making sense out of products of distributions, and integrals of such products, which 
(as argued above) are ubiquitous in quantum field theory, although their exact meaning is 
not usually discussed in standard introductory textbooks. In order to simplify quantum 
field theory calculations, one defines the Wick transform of a power $(/)" := ($", /) so 
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as to satisfy an orthogonality property with respect to Gaussian integration. Recalhng 
the orthogonahty properties of Hermite polynomials (appendix B) and the definition (|^) 
a simple idea is to define the Wick transform in terms of Hermite polynomials: 



:Hfr-c = C{f,fr/'Re^(^-^=j_ 



(17) 



Notice that this depends on the covariance C, and that there is no analogous dependence 
in the analogous harmonic-oscillator definition ([^. 

The orthogonality of two Wick-transformed polynomials is then expressed by 

I rf/ic($) : HJT-c : Hgr-.c = 5^,n n\ mf) $(^?))^J". (18) 

An entertaining exercise is to show this, which we will do in section 3.1.3 (paragraph 2). 
The Wick transform can also be defined recursively by the following equations H] 



Wick transform defined recursively 





■■Hfr-c 


= 1 
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Hfr-.c 


= n:<l>(/)"-i:c 


n = 1,2,... 




Hfr--c 


= 


n = 1,2,... 



where the functional derivative with respect to a distribution is simply 

6 



5$ 



(19) 



Let us check that the Wick transform : $(/)"':(7 defined by eq. ( |19D is the same as the 
Wick transform given in terms of Hermite polynomials in eq. (p^Tj). This is of course to be 
expected, since Hermite polynomials themselves satisfy similar recursion relations, but it 
is a useful exercise to check that it works. To begin with, we establish a property of Wick 
exponentials. 



Wick exponentials 



Let : exp(a$(/)):c be the formal series 



: exp(a$(/)):c = 1 + a: $(/):c + W: HfY:c + 
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where normal-ordering is defined by eq. (|T9|). 
Exercise: Show that 

■ -P'"*'^')- - (exp(U(/)il, ■ (2°) 

One can evaluate the right-hand side of the equation by expanding numerator and de- 
nominator into a power series and dividing one power series by the other. 



oo 



h, t'' 1 °° 



where c„ -f- ^ X]fc=i '^n-fc '^fc ~ = 0. Comparing the resulting series, term by term, to 



the power series expansion of the left side proves equation (|20D . 



Equivalence of Hermite polynomial and recursive definitions 



Exercise: Show the equivalence of (|T^) and ([T7|) . 

Solution: We can explicitly calculate the denominator in eq. (pO[): 

(exp(a$(/)))^^ = j di^ci^) exp{a^{f)) 

n 

= E —/ ^/^^(*) byeq. (0) 

= exp(|a2c(/,/)) byeq. (0). 

Thus, from eq. (^) we find 

: exp(a<|.(/)):c = exp («$(/) - i a'C{f, /)) . (21) 

Multiplying the power series expansions^ of exp(a$(/)) and exp(ia^C(/, /)) and com- 
paring the result term by term to the series expansion of the left side of eq. ( PT| ) yields 

[f] , 

m=0 ■ ^ '' 



•^Similarly to the division of power series mentioned above, the niultiphcation of power series is simply 

/ oo \ / oo \ oo 

^afexM [Y^huxA 



where d„ = YZi=o o-rn bn- 
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Rewriting this expression as 

and using the formula for the defining series of the Hermite polynomials given in appendix 
B, one recovers eq. (p!7|). 

3.1.3 Wick transforms, properties 

Many properties of Wick ordered polynomials can be conveniently derived using the formal 
exponential series. For simplicity, we assume that all physical quantities that one may 
wish to compute (scattering cross sections, etc.) are written with normalization factors of 
1/C{f,g), which in effect lets us set the coincendent-point covariance to unity: C(/, /) = 
1. It can be restored by comparison with eq. (p3|). The properties we are interested in 
are useful exercises to show: 

1. Show that 

: exp ($(/) + $(^)) -.c = exp (-($(/) $(^) ■ exp($(/)):c : exp($(^):c 
Solution: 

: exp($(/)):c7: exp(<l'(^):c 

= exp ($(/) + ^g)) exp (-1 [( ^ff + ( ^{gf 
: exp ($(/) + ^{g)) -.c exp (( $(/) ^{g) , 

where we have used eq. (^Tj) in the first line, and, after completing the square in the 
second factor, again in the second line. Dividing both sides of the equation by the 
second factor completes the demonstration. 

2. Show that 

( : ^{fT-.c : $(^7)'":c ) ,c = ^nm n\ ( $(/) <t>{g) 
Solution: 

If we take the expectation value of both sides of the last line in the proof above, we 
find 

(: exp(<|.(/)):c : exp(<|.((7)):c)^e = ^xp (( $(/) ^g) 

using eq. (p^Q]). Expanding the exponentials on both sides and comparing term by 
term completes the proof. 
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3. Show that 

: = n: ^{fT'^-c - $(/) : ^UT-c 

Solution: 

This is a consequence of the equivalence of Wick ordered functions and Hermite 
polynomials. The expression follows from the recursion relation for Hermite poly- 
nomials given in appendix B. 

4. The definition of Wick transforms given above can be generalized to several fields in 
a very straightforward manner. We quote here some results without proof (details 
can be found in e.g. the book by Simon Q). The reader may find it interesting to 
check that it works: 

: $(A) . . . $(/„+i): = : <l>(/i) . . . $(/„): <l>(/„+i) 

n 

- E /n+l) ^Ul) ■ ■ ■ Hfk-lMfk+l) ■ ■ ■ Hfn): 

k=l 



dfici^):^f,)...^Uy. = 
dficm : $(/i) . . . $(A): : $(^i) . . . ^gj: = forn m 

These latter multi-field expressions reproduce, within this functional framework, what 
is usually referred to as "Wick's theorem" in the creation/annihilation-operator formalism. 
In that formalism, it takes some effort to show this theorem; here we find it somewhat 
easier, relying on familiar properties of the Hermite polynomials. 

3.2 Wick transforms and functional Laplacians 
3.2.1 Definition 

We can also define Wick transforms of functions by the following exponential operator 
expression, which is convenient in many cases (e.g. in two-dimensional quantum field 
theory settings, such as in references |^ and fT^): 

' (24) 



where the functional Laplacian is defined by 



Ac=/d'.dVC(...')^^. (25) 



which, again, depends on the covariance C. 
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Instead of proving eq. (|2^, which is straightforward, we will just illustrate the equiv- 
alence of definition and definition (|T^) in the common example of a 0^ power, when 



the definition reads: 

■■Ayy-c = exp (-1 Ay) ■ 

Expanding the exponential we find 

All higher terms in the expansion are zero. We can now evaluate each term separately; 
the second term is the integral 

d'^x d^x' C(x, x) ^ ^ (p^{y) 
o(p[x) o(p[x ) 

d'^xd'^x' C(x, x)——A(ty^{y) Six' - y) 

0(p[X) 

d'xC{x,y) 12(P\y)6{x-y) 

if the covariance is normalized to unity. We use this result in the third term: 

d'x rfV C{x, x')-^-^ I2ct?{y) 
d(p[x) d(p[T) 

d^x rfV C(x, x')—— 2A(t)(y)6(x' - y) 

0(p{X) 

= 24y" d^xC{x,y)S{x-y) 
= 24 

with the same normalization of the covariance. We collect these results with the appro- 
priate coefficients from the expansion: 

= 0^(l/)-^-12 0%) + ii-24 

= ^\y)-6^^(y) + 3 

= He4(0(i/)), 
which completes the illustrative example. 
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3.3 Further reading 



Although we hope to have given some flavor of some of the techniques and ideas of quan- 
tum field theory mathematical-physics style, we have of course really only given a few 
illustrative examples and demonstrated some simple identities. For more on the math- 
ematical connection between Wick transforms on function spaces and Wick ordering of 
annihilation and creation operators, we recommend textbooks such as and ||^. 



3.4 An application: specific heat 

In this last section, we discuss one example of a physical application of some of the 
results above. By using the connection (0) between the Wick transform and Hermite 
polynomials, we show how one can exploit standard properties of those polynomials to 
simplify certain calculations. 

Consider the familiar generating function of Hermite polynomials (but for the scaled 
polynomials (|)): 

e""-"'/2 = ^He„(a;) aVn! . (26) 

This generating function gives a shortcut to computing some quantum effects in two- 
dimensional quantum field theory, where normal orderering is often the only form of 
renormalization necessary. An important issue is the scaling dimension of the normal- 
ordered exponential : e*^*^: , where p is a momentum. (The real part of this operator can 
represent the energy of a system where is the quantum field.). In other words, the 
question is: if we rescale our momentum p Ap, equivalent to a rescaling x —>■ A~^x in 
coordinate space, how does this operator scale? Since an exponential is normally expected 
to be dimensionless, one might guess that the answer is that it does not scale at all, i.e. 
that the scaling dimension is zero. In fact, this is not so due to quantum effects induced 
by the normal ordering. Passing to the functional integral, we can easily compute the 
effect of the normal ordering (here, the Wick transform): 

oo ^ 

.^ip4>,^ = y^-Atpr:r--c 

n=0 
oo ^ 

= E;^(^P)"^"^'He„(0/v^) 

n=0 

= e^p'^e^P"^, (27) 

where we used the definition (|T^ and the previously given generating function (pB]). Now, 
the covariance (Green's function) C is a logarithm in two dimensions, i.e. the solution 
of the two-dimensional Laplace equation is a logarithm. To regulate divergences when 
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p — > oo, one introduces a cutoff A on the momentum, which makes C = In A. Substituting 
this into ( P?] ) yields the answer 

: e'^'^: = AP'/^e'^'^ . (28) 
Thus, the anomalous scaling dimension, usually denoted by 7, is 7 = for the expo- 



nential operator. This is an important basic result in conformal field theory (see e.g. |TT 
p. 451). Here, the just comes from the a^/2 in the generating function (126|)! 

How could such a quantum effect be measured? Consider the "two-dimensional Ising 
model with random bonds" (p. 719, [l^). This is just the familiar Ising model, but 
one allows the coupling between spins to fluctuate, i.e. the coupling becomes a space- 
dependent Euclidean field in the spirit of previous sections. The energy of the system 
is described by (the real part of) an exponential operator : e*^*^: as stated above. Using 
renormalization group methods, the anomalous dimension (^8]) leads to a formula for the 
specific heat in this system (eq. (356) in |T^, where the derivation is also given). The 



specific heat is in principle directly measurable as a function of the temperature, or more 
conveniently, as a function of ^ = (T— Tc)/Tc, the dimensionless deviation from the critical 
temperature. The renormalization group description predicts a certain double logarithm 
dependence on 9 that could not have been found by simple perturbation theory, and it 
uses as input the result (p^). 

Admittedly, the telegraphic description in the previous paragraph does not do justice 
to the full calculation of the specific heat in the two-dimensional Ising model with random 
bonds. Our purpose here was simply to show how the Wick transform reproduces the 
quantum effect (pSD, and then to give some flavor of how this effect is measurable. 



4 Conclusion 

In this article we have shown that the scope of normal ordering has expanded to set- 
tings beyond the original one of ordering operators. Several different definitions of Wick 
ordering of functions have been discussed and their equivalence established. 

For deeper understanding and further applications of these ideas, the interested reader 
is invited to consult the quoted literature, which is a selection of texts we found particu- 
larly useful. Specifically, for the physics of functional integrals we enjoy 0. For a more 
mathematically oriented treatment we find quite useful. 
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A Wick ordering of operators in QFT 

In this appendix, we briefly remind the reader how the need for Wick ordering arises in 
the operator formulation of quantum field theory. All of this is standard material and can 
be found in any introductory book on quantum field theory (e.g.[|]), albeit in lengthier 
and more thorough form. We set ^ = 1 throughout. 

Consider a real scalar field x) of mass m defined at all points of four- dimensional 
Minkowski spacetime and satisfying the Klein-Gordon equation 

g2 



V^ + m^ 0(t,x) = 0. 



The differential operator in the parenthesis is one instance of what we call D in the text. 
The classical Hamiltonian of this scalar field is 

H=IY1 [(^(^' ^))' + ^))' + ^)] (29) 

X 

where vr is the variable canonically conjugate to (p, in fact it is simply vr = d(j)/dt. Here 
we can think of the first term as the kinetic energy, and the second as the shear energy. 
This classical system is quantized in the canonical quantization scheme by treating the 
field (j) as an operator, and imposing equal-time commutation relations 

[0(t,x),0(t,xO] = 0, 
[7r(t,x),7r(t,x')] = 0, 
[0(t,x),7r(t,x')] = «53(x-x'). 

The plane-wave solutions of the Klein-Gordon equation are known as the field modes, 
uy_{t, x). Together with their respective complex conjugates ul_(t, x) they form a complete 
orthonormal basis, so the field (p can be expanded as 

(j){t, x) = ^[flk ui^it, x) + al uKt, x)]. 

k 

The equal time commutation relations for (p and vr are then equivalent to 

[flk , ctk'] = , 



[al,al,] = 



k ; "k' 
Ctk , oj^/] = (^kk' 
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These operators are defined on a Fock space, which is a Hilbert space made of ra-particle 
states {n = 0,1, . . .). The normalized basis ket vectors, denoted by | ), can be constructed 
starting from the vector |0), the vacuum. The vacuum state |0) has the property that it 
is annihilated by all the operators: 

ak|0) = 0, Vk. 



In terms of the frequency u!]^ = CA/|kp + m^, the Hamiltonian operator obtained from 
H = ^ ^(4ak + akal)ujk = 5]](4^k + |) ^^k , 

k k 

where in the last step we used the commutation relations from above. Calculating the 
vacuum energy reveals a potential problem: 



|^|0) = (0|0)5^1a;k = $^i^k oo 



where we have used the normalization condition (0|0) = 1. This infinite constant can be 
removed as described in the text. 

Propagation amplitudes in quantum field theory (and hence scattering cross sections 
and decay constants) are given in terms of expectation values of time-ordered products of 
field operators. These time-ordered products arise in the interaction Hamiltonian of an 
interacting quantum field theory. The goal is to compute propagation amplitudes for these 
interactions using essentially time-dependent perturbation theory, familiar from quantum 
mechanics. At leading order in the coupling constant, these products can be simplified, 
and the zero-point constant energy removed by using Wick's theorem. 

The way Wick ordering is applied in practice to calculations in QFT is through "Wick's 
theorem", which gives a decomposition of time-ordered products of field operators into 
sums of normal-ordered products of field operators (again, we refer to e.g. p). In this 
paper, in section p.l.3| . Wick's theorem appears in the functional-integral formulation of 
the theory. 



B Properties of He^(x) 

Here we list a few useful properties of the scaled Hermite polynomials. More can be found 
ini. 



Defining series: 

[i] 

Rejx) = y"(-ir—, — 7 



n-2m 



m=0 
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where [n/2] is the integer part oin/2. 
Orthogonahty: 



Generating function: 

Recursion relation: 
The first five: 



/oo 
dx e"^^/^ He„(a;) He,„(x) = bmn'^f^ 
-oo 



exp {xa — \ o?^ = B.en{x) — - . 

n=0 



He„+i(x) = xHe„(a;) - nHe„_i(x) . 



n! 



Heo(a;) 


= 1, 




Hei(a;) 


= X , 




He2(x) 


^ x'- 


1, 


He3(x) 


^ x'- 




He4(x) 


^ x'- 


6x^ + 3 
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